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I. INTRODUCTION 



It was speculated by Witten that some of the dark matter in the universe might 
possibly exist in the form of strange quark matter (SQM), which consists approximately of 
the same amount of up, down and strange quarks, and might be formed after the big bang 
when the universe underwent the quark to hadron phase transition. Another place that 
people can expect to find SQM is in the "little bang", i.e. high energy heavy-ion collisions. 
It was supposed that if a hot quark-gluon plasma (QGP) is formed after the collisions, 
the cooling process of the QGP might result in droplets of SQM, so called strangelets, which 
could serve as an unambiguous signature for the QGP formation. In fact, some heavy-ion 
experiments @] are searching or planning to search for strangelets, hence reliable theoretical 
calculations of strangelets formation in heavy-ion collisions are in urgent need. 

Recently Lee and Heinz has presented a comprehensive study of the phase structure 
of a bulk strange quark matter in equilibrium with a bulk hadronic gas. Assuming Gibbs 
phase equilibrium and an isentropic expansion of the hot fireball, they found under some 
conditions strange quark matter can actually survive the hadronization and cooling process, 
giving rise to strangelet formation. On the other hand, taking into account the particle 
evaporation from the surface of the fireball, some rate calculations have been carried 
out to study the distillation and survival of strangelets in heavy-ion collisions. In this work 
we study finite size effects on strangelet formation, an important aspect which has been 
neglected in the above mentioned calculations. 

It has been shown that finite size effects may strongly destabilize small strangelets at 
zero temperature 0, and significantly modify the properties of strangelets at finite tem- 
perature [§]. Finite size should be important as well in the study of the possible formation 
of strangelets, which are expected to be produced small in heavy-ion collisions The 
present paper follows the work of Lee and Heinz to study the phase diagram and evolution 
of strangelets in heavy-ion collisions, while extending it to include the essential finite size 
effects. We find that there are indeed significant finite size modifications of the phase dia- 
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gram and evolution of strangelets. A realistic calculation that attempts to provide reliable 
information for the strangelets searches should therefore take finite size effects into account. 

To study SQM in the context of the MIT bag model we have generally two ways to 
account for finite size effects. One employs the shell model, namely solving the Dirac equa- 
tion with the bag boundary conditions and populating the energy levels of the bag P,|5|,pi? 



However, shell model calculations are very tedious even at zero temperature, and become 
extremely difficult when one tries to study the phase diagram and evolution. Another way to 
introduce finite size effects of strangelets is to explore a continuous density of states within 
the framework of multiple reflection expansion which leads to the liquid drop model 
||7|JT^. As shown by Madsen the liquid drop model reproduces the overall structure of 
shell model results for strangelets, and allows much easier calculations than shell model. In 
this work we shall use the liquid drop model to incorporate finite size effects into the studies 
of phase diagram and evolution of strangelets. 

This paper is organized as follows. In Sect, y we give a description of equations of state 
for the QGP and hadronic phase, with an emphasis on the inclusion of finite size effects in 
the QGP phase. In Sect. |T| we study the phase diagram of a strangelet in equilibrium with a 
finite hadronic phase, assuming the conservation of total baryon number and strangeness in 
the system. It is shown that the inclusion of finite size effects can modify the phase diagram 
in a significant way. The isentropic expansion of the system and possible formation of cold 
strangelets are explored in Sect. We find that in contrast with bulk case the inclusion of 
finite size effects will give even more restrictive conditions under which dynamical formation 
of a cold strangelet is possible. We summarize our work in Sect. 0. 



II. EQUATION OF STATE 



3 



A. The QGP Phase 



We consider the QGP phase as a gas of free massless up and down quarks (m^ = = 0), 
massive strange quarks (m^ = 150 MeV), their antiquarks and gluons confined in an MIT 
bag. In the multiple reflection expansion of Balian and Bloch [n| the density of states for 
particle species i present in the bag is given by, 

For spherical QGP droplets = 47r(i?'5)^/3 is the volume of the bag, S"^ = 47r(i?*5)2 jg ^.j^g 
surface area, and = 8ttR^ is the extrinsic curvature of the bag surface. The factor Qi is 
the statistical weight (6 for quarks and antiquarks, and 16 for gluons). 

The functions fg^ and f^^ are surface and curvature terms for particle species i. The 



surface term for quarks was given by Berger and Jaffe [111] as 



In particular, the surface term for massless quarks and gluons is zero. The curvature term 
for gluons is f^f^ = — l/Gvr^ |l^, while for massless quarks fc\^) = — l/247r^ ||15|. It has 



been shown by Madsen [0] that the following ansatz works for the curvature term for massive 
quarks: 

The thermodynamic potential of the QGP phase can be derived from, 

Q = Y.Q, + BV'^, (4) 

i 

where B is the bag constant, and the thermodynamic potential of particle species i can be 
written as 

Q, = TT dk—^ ln[l ± exp{-{E,{k) - fi,)/T)] 
Jo aft 

= n.yV'^ + n.^sS'^ + fi,,cC«, (5) 



with the upper sign for fermions and the lower for bosons. Here /i^ is the chemical potential. 
We have assumed isospin symmetry so that yU„ = /i^ = fiq, and there are two independent 
chemical potentials in the QGP phase, namely fiq and /ig. To get an impression of the 
finite size effects introduced here, we have plotted in Fig. 1 the temperature dependence 
of individual surface and curvature contributions of particle species to the thermodynamic 
potential. One sees from Fig. 1 that they are all increasing functions of the temperature. 
For a baryon and strangeness free QGP droplet {i.e. fiq = fig = 0) the surface and curvature 
contributions vanish at zero temperature (Fig. 1(a) and (b)), while for nonvanishing chemical 
potentials there are finite contributions even at zero temperature. We see that significant 
finite size effects can be given by an approach of multiple reflection expansion, which can be 
relevant in the study of QGP phase structure. 

After the construction of the thermodynamic potential, we can readily obtain the ther- 
modynamic quantities of the system as follows. The number of quarks can be derived from 

N? = - (f^) , (6) 



and 



N? = \Y.N^ (7) 



gives the total baryon number of the QGP droplet. The total entropy of the QGP phase is 

'dVt\ 
df I 



B. The Hadronic Phase 



Now we consider the equation of state for the hadronic phase. We model the hadronic 
phase as a mixture of Bose-Einstein and Fermi-Dirac gases of mesons vr, K, t], baryons 
A^, A, S, S, Q, A and their anti-particles p!6| , p!2| . We will show below that using this 



limited particle spectrum gives practically identical results compared with the case in which 
a spectrum of all hadrons up to M ~ 2 GeV is used, thus sufficient for the study of the phase 



diagram. We can write down the pressure, particle number, and energy density expressions 
for pointlike hadron species i as follows, 

1 
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671^ Jo Ei(fc)e(^«('=)-^>)/^±l 

2712 70 " e(^«W-'^>)/^±l' ^ ^ 

where gi =(spin) x (isospin) is the degeneracy factor for each hadron, the ± signs correspond 
to Fermi-Dirac and Bose-Eistein statistics, /ij is the chemical potential of hadron species i 
and given by the net numbers of light quark and strange quark, AiV*, AA^*, as 

/X. = AiVj /xf + ANl /xf , (12) 

if chemical equilibriums are assumed. 

The thermodynamic quantities for the hadronic gas are corrected by a Hagedorn eigen- 
volume factor, (1 + e'^^ /4B)~^, to take into account the repulsive interactions between the 
hadrons [17,1^, which is essential for the existence of phase transition in the considered 



model, 

= Kt— y biNf, (14) 

= "E s^Nf\ (15) 

1 + £pV45Y 

where = Z^j^f*, and 6j, Si are the baryon number, strange valence quark number of 
hadron species i respectively. 



We want to emphasize that a finite hadronic phase, instead of a bulk one |T^, is con- 
sidered in this work. In the case of a bulk QGP phase in equilibrium with a bulk hadronic 
phase, the conservation of total net baryon number and strangeness in the system is ex- 
pressed in terms of a fixed strangeness fraction, namely the net number of strange quarks 
per baryon. In our work, however, to include the finite size effects we have considered a finite 
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QGP phase, whose volume, baryon number and strangeness have been specified. Therefore, 
if the system consisting of a QGP phase and a hadronic phase is assumed to conserve total 
net baryon number and strangeness as a whole (the conservation of total net baryon number 
and strangeness will be discussed furthermore below), the volume of the hadronic phase has 
to be specified as well in order to conserve the total net baryon number and strangeness in 
a consistent way. On the other hand, we have not included surface and curvature contribu- 
tions from hadrons in our calculations since we have not found appropriate expressions for 
them |T^. Nevertheless, there were expections that these contributions from hadrons may 



be less important compared with those from quarks and gluons 1 19,20 



III. PHASE DIAGRAM WITH FINITE SIZE EFFECTS 

We start this section with a description of our scenario. If a baryon-rich QGP is formed 
in high energy heavy-ion collisions with an equal number of s and s quarks, the system 
will suffer a strangeness enrichment due to the early black-body radiation of more kaons 
(qs) than anti-kaons (qs) off the fireball, and then the system reaches the phase coexistence 
region We shall study the phase diagram (this section) and further phase evolution (next 
section) of the resulting system with finite net strangeness. 

Following the work of Lee and Heinz we consider a quasistationary phase transition 
between the QGP and hadronic gas, which assumes that the thermodynamic equilibrium 
time as well as the hadronization time are small compared to the lifetime of the mixed 
phase so that a complete thermodynamic equilibrium can be held during the whole phase 
transition. In this work we have not taken into account the effects of particle evaporation 
off the system which may change the total net baryon number, strangeness content of the 
system 0,0. As a result, the total net baryon number in the system is kept a constant. 
Furthermore, since we expect that the typical time scales in heavy-ion collisions are too 
short to establish a flavor equilibrium due to the weak interactions which can convert strange 
quarks to nonstrange quarks and vice versa, it is assumed that the total net strangeness in 
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the system is conserved. Note that the total net baryon number and strangeness can be 
simultaneously conserved in a consistent way only when a finite hadronic phase is considered 
in connection with a finite QGP phase, as we discussed in Sect. 

The Gibbs equilibrium between the QGP and hadronic phase is given by pQ = 
(mechanical equilibrium), T'^ = (thermal equilibrium), /i^ = /i^, fif = fi^ (chemical 
equilibrium). Taking generally the temperatures and chemical potentials in the two phases 
equal, we are left with the pressure balance between the two phases, 

P'^{T,fi„fi,-R'^) = P^(T,/i„/i,). (16) 

Note that the pressure generated by the QGP phase is a function of the QGP droplet radius 
since we have included finite size effects into the equation of state for the QGP phase. 
The conservation of total net baryon number and strangeness in the system gives two 
constraints, 

and 

TVf = iVf (T, /i„ /i,; R^) + iVf (T, /i„ /i,; V""), 
where the volume of the hadronic phase is given by 

yH ^ ytot _ yQ_ ^-^g^ 

The phase diagram of the system is determined by Eqs. ([T6| ) — ([T9| ) , and shown in Fig. 2(a) 
for parameters i?^/^ = 180 MeV, = 150 MeV, and a fixed strangeness fraction (namely net 
number of strange quarks per baryon) = 0.5. Full lines in Fig. 2(a) illustrate our results 
for a finite system with total net baryon number A*"^"* = 100 and strange quark number 
A^*°* = 50. To be more specific, the full line separating the pure QGP phase and mixed 
phase (MP) is obtained by setting = in Eqs. ([T^) and (0), and solving Eqs. (|IBp -([!§) 
for (/iq, R^) at a given T, A*"^"* and A^*°*. On the other hand, drawing the boundary 
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(18) 



between the mixed phase (MP) and hadronic gas (HG) is a httle more dehcate. For a given 
value of we solve Eqs. (|1|)-(|1|) for (/i,, /i„ V^) at fixed T, A^*"* and A^f*. We 
then increase the value of until is so large that Eqs. (0) — (^) have no solutions 
for [fig, fis, 1 V^), which implies that the QGP and hadronic phase can no longer coexist 
in equilibrium, and thus we have reached the boundary between the mixed phase (MP) 
and hadronic gas (HG). We observe that at the MP— HG boundary the QGP droplet may 
have non- vanishing baryon number and strangeness, which is of course a consequence of the 
inclusion of finite size effects. To check the consistency of this solution procedure, we have 
also tried to decrease the value of V"*°* while solving Eqs. ([161) — (|T9|) until no solutions exist, 
and found that we exactly arrived at the QGP— MP boundary, which we have obtained 
above. 

In comparison with the above case of a finite system, the dashed curves in Fig. 2(a) show 
the phase diagram of a bulk QGP in equilibrium with a bulk hadronic phase, i.e. without 
finite size effects, which is identical to the results of Lee and Heinz in Fig. 6(c) of Ref. One 
sees from Fig. 2(a) that the hadronic phase grows significantly with the inclusion of finite 
size effects, while the QGP— MP boundary is little affected by the finite size corrections. 
We shall show in next section that this growth of the hadronic phase will greatly reduce the 
survival probabilities of cold strangelets which we wish to find in heavy-ion collisions. 

The dotted lines in Fig. 2(a) is plotted for A^"^"* = 50, A^*°* = 25, to show how the finite 
size modifications of phase diagram vary with the given iV^"* and A^*°*, in comparison with 
the bulk case (dashed lines) and the case of Nl°^ = 100 and A^*°* = 50 (full lines). It is of no 
surprise that with smaller iV^°* and A^*°* the smaller system has its phase diagram modified 
more prominently by the finite size effects. 

To see how the modifications of phase diagram due to the inclusion of finite size effects 
are affected by the bag constant B, we have plotted in Fig. 2(b) and 2(c) the phase diagrams 
for i?^/^ = 145 MeV and 235 MeV, fs = 0.5, with all dashed lines for a bulk system, and 
all full lines for a finite system with total net baryon number A*"^"* = 100, strange quark 
number A^*°* = 50. We see that finite size modifications are more significant for larger bag 



constants. This can be understood roughly that for given A*"^"*, A^*"*, and T, larger bag 
constants generally correspond to more compact and smaller QGP droplets in equilibrium 
with a hadronic gas, and thus more important finite size effects. However, we want to point 
out that the full line in Fig. 2(c) separating the mixed phase (MP) and hadronic gas (HG) 
is somehow out of control at regions of high baryon number density pf,. This difficulty 
may come form two origins. First, at regions of high baryon number density the chemical 
potentials may become very large and lead to some numerical uncertainty in the relevant 
Bose-Einstein integrations, which can be solved by further careful examinations. Another 
possibility is that at these regions the radius R'^ of QGP droplet in equilibrium with a finite 
hadronic gas tends to become so small (^ 1 fm) that the multiple refiection expansion in 
powers of which we have adopted to introduce finite size effects, could be somewhat 

dangerous, as previously warned by Madsen |^. Fortunately, in our calculations this kind of 
disease sets in generally for very large bag constants. As it is expected that for sufficiently 
large bag constant strangelets can neither be stable nor metastable even at zero temperature 
and the formation of cold strangelets in heavy-ion collisions is very unlikely [0,0, therefore 
this problem will not become serious if we are only interested in the cases in which the 
formation of a cold strangelet can be possible. 



IV. FORMATION OF COLD STRANGELETS 

In this section we turn to study the evolution of the system to see how a cold strangelet 
could be formed in the final stage of evolution and what its characteristics look like. While 
a full dynamical investigation of the fireball evolution is beyond the scope of this work, we 
assume that the system expands adiabatically, conserving its total entropy during the whole 
evolution p6|pl| . Note again that we have not included the process of surface evaporation 



which may change the entropy, net baryon number and strangeness of the system. Therefore 
what we are exploring is a system that experiences a smooth hydrodynamic expansion, 
conserving the total entropy S**"*, net baryon number Nl°^, and strangeness A'^^"* [Bl. As we 
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have already described in Sect. [IT4 the system considered here may have nonvanishing net 
strangeness due to the early surface radiation of hadrons after the collisions. 

We first look at a QGP fireball which is expanding isentropically and cooling down until 
it reaches the boundary of mixed phase. The isentropic expansion trajectories in the pure 
QGP phase are determined by solving 

= iV?(T,/i„/i,;i?'3), (20) 
iVf* = iV«(T,/i„/i,;i?«), (21) 

and 

for {T, fig, fig) at fixed Nj^"^, A^*°*, 5**°* and various R^. Increasing the value of until the 
pressure balance Eq. (0) is possible, we will arrive at the boundary between the QGP and 
mixed phase, which we have obtained in the phase diagram given in last section. 

Afterwards the isentropic expansion trajectories go into the mixed phase and are given 
by Eqs. (0) — (|T^, together with 

gtot ^ 5Q(T,;i„^,;i?Q) + S^(T,/x„/i,;\/^), (23) 

which are solved for (T, fig, fig, R^ ^V^) at fixed A*"^"*, A"*"*, 5**°* and various V^°^. When we 
increase the value of V^°^ while solving these equations, two phenomena may be observed. 
One is that these equations have no solutions for [T, fig, fig, R^, V") when V^"^ is large 
enough, which is found to be the very case in which the expansion trajectory hits the 
MP— HG boundary at some point. This observation further confirms the consistency of our 
treatment of finite size effects. In this case the expansion trajectory will then enter the 
hadronic phase, corresponding to a complete hadronization, and no cold strangelet could 
be formed. Another important phenomenon is that for V^°* 00 Eqs. (|T6D — (|T9|) and ( |23D 
can always be solved at fixed A"^"*, A^^"* and 5**°*, giving T ^ 0. This implies that a cold 
strangelet could appear in the final stage of the evolution, which is of great interest in 
heavy-ion experiments. 
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As manifestations of the above discussions, we show in Fig. 3 isentropic expansion tra- 
jectories of the system through the phase diagrams for B^^^ = 180 MeV, A^"^"* = 100, and 
Nl°* = 0, 50, 100, 200, or strangeness fractions fs = 0, 0.5, 1.0, 2.0. Our Fig. 3 can be 
compared with the results of Lee and Heinz for bulk phases given in Fig. 6 of Ref. ^ to see 
how the inclusion of finite size effects can modify the isentropic expansion trajectories of the 
system. It can be seen in Fig. 3 that for most cases the expanding system will hadronize 
completely. In particular, for fs = 0.5 and S^°*/Nl°^ = 5 (Fig. 3(b)), in contrast with the 
bulk case (Fig. 6(c) in Ref. P|) in which the system will always stay inside the mixed phase 
and give rise to cold strangelet formation, one finds that the inclusion of finite size effects 
tends to make the formation of cold strangelets even more difficult. In fact, with an impres- 
sive observation from Fig. 2(a) on how the phase diagram can be drastically modified by 
the finite size of the system, we can expect that currently existing calculations of possible 
formation of cold strangelets in heavy-ion collisions without finite size effects may have to 
become somehow less optimistic with the inclusion of finite size effects. 

The inclusion of finite size effects may also affect in some way the strangeness separation 
process during the phase transition which is an essential ingredient for the possible formation 
of cold strangelets in heavy- ion collisions Fig. 4 shows the evolution details of the QGP 
droplet corresponding to the isentropic expansion trajectory for a system with zero net 
strangeness and S'*°*/^b°* = 5 in Fig. 3(a). In Fig. 4(a) the system experiences a cooling in 
the pure QGP phase, and reheating in the mixed phase which is mainly due to the energy 
surplus when quarks in the QGP are converted into hadrons in the hadronic gas. Fig. 4(b) 
shows a smooth variation of chemical potentials fig and fig during the phase transition. In 
particular, /i<j remains zero in the pure QGP phase whose net strangeness is zero, but increase 
continuously in the mixed phase even though the net strangeness of the system remains zero, 
indicating the occurrence of strangeness separation which can be seen also in Fig. 4(c) from a 
continuous increase of the strangeness fraction of the QGP droplet. During the evolution 
in the mixed phase the net baryon number A*"^ (Fig. 4(d)) and radius (Fig. 4(e)) of the 
QGP droplet decrease monotonously. However, at the boundary between the mixed phase 
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and hadronic gas A*"^ and R'^ do not vanish, and only a rather moderate value of can 
be reached. This is quite different from the corresponding cases for bulk phases studied 



previously pj2^, in which the QGP phase diminishes at the MP— HG phase boundary and 
the strangeness separation mechanism works in a much more effective way, leading to a 
much larger value of f^. In fact, as we have already known in the studies of phase diagrams 
that the hadronic phase grows on account of the mixed phase with the inclusion of the 
finite size effects, we can expect that the inclusion of the finite size effects tends to make the 
strangeness separation proceed to a less extent. Since the strangeness separation mechanism 
appears so important in calculations of the possible formation of cold strangelets in heavy- 
ion collisions, we hereby feel it a substantial task to include finite size effects in the future 
relevant investigations. 

Now we study the formation of cold strangelets in heavy-ion collisions. The isentropic 
expansion trajectories in Fig. 3(c) and (d) indicate that for S*°^ /N^"^ = 5 the system will 
always stay in the mixed phase and expand infinitely to zero temperature and zero baryon 
number density, resulting in the formation of a cold strangelet. As an essential consequence 
of the inclusion of finite size effects, we can now obtain the properties of the surviving 
strangelet. For example, we list in Table. | the properties of the resulting cold strangelet 
corresponding to the expansion trajectories for S^"^ /Nl°*'=5 and /*°* = 1.0 (Fig. 3(c)), 2.0 
(Fig. 3(d)) respectively. Note in the case of /*°* = 1.0 all net strangeness is in the droplet 
^jyQ _ 100), which is due to the fact that at T ^ chemical potentials fig, fig in the phase 
equilibrium configuration acquire values such that only nucleons but no strange baryons are 
present in the hadronic gas. 

However, the properties of cold strangelets given in Table. | should not be taken too 
seriously for two reasons. First, an isentropic and equilibrium expansion of the fireball all 
the way to zero temperature assumed in this work is not a realistic scenario. Effects of 
particle evaporations off the surface of the fireball can be very important for the strangelet 
formation MM. At sufficiently low temperature and density the QGP droplet may decouple 
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from the surrounding hadronic gas and the thermodynamic equihbriums between them may 
be broken. On the other hand, the hsted properties of the strangelets rely heavily on 
the parameters such as bag constant B, strange quark mass rris, and on the given initial 
conditions of the system in terms of A*"^"*, A^*°* and 5**°*. Nevertheless, Table. | can give us 
a first impression of the characteristics of cold strangelets resulted from a simple dynamical 
evolution, and may provide clues for future construction of realistic dynamical evolution 
models. 

The isentropic expansion trajectories in Fig. 3 imply that only for sufficiently low 
S^°^ /Nl°^ or large fs does the system allow the formation of cold strangelets. Following 
Lee and Heinz we have plotted in Fig. 5 the maximum total entropy per baryon S^°^ /Nl°* 
permitted by the formation of cold strangelets versus the strangeness fraction fg, for two bag 
constants B^/^ = 145 and 180 MeV, and a finite system with net baryon number N^"^ = 100. 
In comparison with the results for bulk phases in Fig. 8 of Ref. 0, we see that the inclu- 
sion of finite size effects give much more restrictive conditions under which cold strangelets 
could be formed. Even for bag constant as small as B^^'^ = 145 MeV, which is expected to 
be the minimum value compatible with the usual nuclear matter phenomenology, and for 
strangeness fraction /<, = 3.0, an upper limit S^°^ /Nl°^ ^ 52 is required to make the for- 
mation of cold strangelets possible. However, these results may depend on the net baryon 
number A^^°*, or the size of the system, and are to be made more convincing by further 
detailed investigations. 

V. SUMMARY 

In the framework of multiple refiection expansion we have studied finite size effects on 
the phase diagram and evolution of a strangelet in equilibrium with a finite hadronic gas. 
Our aim is to check how important the finite size effects can be in the studies of possible 
formation of cold strangelets in high energy heavy-ion collisions. We have shown that there 
are indeed significant finite size modifications of the phase diagram, which are found to 
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become more important for a system with larger bag constant B or smaller total net baryon 
number N^"^. 

Assuming an isentropic expansion of the system we have studied the finite size effects 
on the evolution of a strangelet. It is observed that the inclusion of finite size effects tends 
to make the strangeness separation mechanism less effective, and will give very restrictive 
conditions for the cold strangelet formation in heavy-ion collisions, expressed in terms of 
windows of the bag constant B, the total entropy per baryon S^°* /Nl°^, and the strangeness 
fraction /*"*. With the inclusion of finite size effects we have been able to obtain quantita- 
tively the properties of cold strangelets emerging in the final stage of the isentropic expansion 
of the system. 

Even though the treatment of finite size effects within the framework of multiple reflection 
expansion may fail for very small strangelets as we have uncovered in the studies of the phase 
diagrams, and the isentropic and equilibrium evolution scenario studied in this work is not 
a realistic one, yet we can still conclude after this work that finite size effects do play 
an essential role and should be included in future realistic investigations of the possible 
formation of cold strangelets in high energy heavy-ion collisions. 
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TABLES 

TABLE I. Properties of surviving strangelets for a system with B^/^ = 180 MeV, 

nis = 150 MeV, A^"* = 100, and S*°^/Nl°^ = 5. The second column describes a cold strangelet 
resulting from the isentropic expansion trajectory S^°^/N^°^ = 5 in Fig. 3(c) with A*°* = 100 or 
equivalently = 1.0, and the third column corresponds to the isentropic expansion trajectory 
Siotjj^tot = 5 in Fig. 3(d) with Af* = 200 or /*°* = 2.0. 

* = 1-0 /f * = 2.0 

(fm) 3.07 3.43 

A^^ 62.7 84.2 

= N^/V^ (fm-3) 0.52 0.50 

100. 170. 

fQ = NQ/N^ 1.59 2.02 

Z/N^ = (1 - /F)/2 -0.30 -0.51 

E'^/Nj^ (MeV) 1184 1250 
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FIGURES 

FIG. 1. Temperature dependences of the surface and curvature contributions to the thermody- 
namic potential of ligiit quarks ("q"), massive strange quark ("s"), and giuon (''g'') for vanishing 
and nonvanishing chemical potentials. The "sum" in (b) and (d) stands for the total curvature 
contributions. 

FIG. 2. (a) Phase diagrams for a QGP phase in equilibrium with a hadronic gas (HG), with 
bag constant S^/^ = 180 MeV and fixed strangeness fraction fs = 0.5: dashed lines for a bulk QGP 
and a bulk hadronic gas; full lines for a finite system with total net baryon number A^^°* = 100 
and net strange quark number AT*"* = 50; dotted Hues for N^"^ = 50 and iV*"* = 25. Note that 
the pure QGP phase and the hadronic gas (HG) are separated by a finite mixed phase (MP), (b) 
Phase diagrams for a QGP phase in equilibrium with a hadronic gas (HG), with B^/^ = 145 MeV 
and fixed fs = 0.5: dashed lines for a bulk QGP and a bulk hadronic gas; full lines for a finite 
system with AT*"* = 100 and ATf* = 50. (c) The same as in (b) but for = 235 MeV. 

FIG. 3. Isentropic expansion trajectories of a finite system with various total entropy per baryon 
S^°^/Nl°^, total net baryon number N^"^ = 100, and net strange quark number (a) A^*°* = 0; (b) 
ATf* = 50; (c) ATf * = 100; (d) iVf * = 200. (SV^ = igo MeV, = 150 MeV) 

FIG. 4. Variations of (a) temperature T, (b) chemical potentials /Xg, fig, (c) strangeness fraction 
/P, (d) net baryon number , and (e) radius R'^ of the QGP droplet during the isentropic 
expansion of a system with S^°^ /N^°^ = 5, N^"^ = 100 and A^*'** = (corresponding to the isentropic 
expansion trajectory S^°*/Nl°* = 5 in Fig. 3(a)). (5^4 = igO MeV, rus = 150 MeV) 

FIG. 5. Maximum total entropy per baryon S^°^ /N^"^ permitted by the formation of cold 
strangelets as a function of the total strangeness fraction /*°* for a system with Nl°^ = 100 and 
two values of bag constant B. 
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